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Abstract. - We propose a simple surface potential describing the planar degenerate anchoring 
of nematic liquid crystals, i.e., the tendency of the molecules to align parallel to one another 
along any direction parallel to the surface. We show that, at lowest order in the tensorial 
Landau-de Gennes order-parameter, fourth-order terms must be included. We analyze the 
anchoring and wetting properties of this surface potential. In the nematic phase, we find the 
desired degenerate planar anchoring, with positive scalar order-parameter and some surface 
biaxiality. In the isotropic phase, we find, in agreement with experiments, that the wetting 
layer may exhibit a uniaxial ordering with negative scalar order-parameter. For large enough 
anchoring strength, this negative ordering transits towards the planar degenerate state. 



Introduction. - Nematic liquid crystals are non-polar liquids consisting of elongated 
molecules aligned, on the average, parallel to a common direction ±n [1]. The unit vector n 
is called the director. Exceptionally, in the case of "negative order," the molecules may be 
distributed perpendicularly to n. The nematic order is quantified by a symmetric traceless 
tensor Q, of Cartesian components Qij {i,j — 1,2,3). Microscopically, calling m a unit 
vector parallel to the long axis of the molecules, the tensor order-parameter can be defined as 
Qij = {rriimj) — ^Sij, where the brackets denote statistical averaging and Sij is the Kronecker 
symbol. For uniaxial ordering, which corresponds to the preferred state of the nematic in the 
bulk, two of the eigenvalues of Q are equal; in this case Qij can be written as S{ninj — ^Sij), 
where — ^ < S' < 1 is the scalar order-parameter. Note that S = corresponds to the isotropic 
phase (randomly oriented molecules), occurring at temperatures above the nematic-isotropic 
phase transition; 5* > corresponds to the ordinary nematic, in which the molecules are 
aligned on the average parallel to n, while S < corresponds to the "negative" case, in which 
the molecules are on the average orthogonal to n. 

The alignment of a nematic liquid crystal by a bounding interface is of considerable in- 
terest both for fundamental and technological reasons. In the most common situations, this 
"anchoring" effect is well described by a surface free energy depending only on the director 
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n at the interface [1]. Three typical situations are: i) homeotropic anchoring, where the 
preferred, or "easy" , average orientation corresponds to n normal to the interface, ii) planar 
anchoring, where the preferred average orientation corresponds to n lying in one particular 
direction parallel to the interface, and in) planar degenerate anchoring, where all the planar 
orientations for n are equivalent easy directions. The latter situation is often realized with 
isotropic bounding fluids, e.g., glycerin [2] or gallium [3]; degenerate anchorings can also be 
obtained with polymer-grafted surfaces [4]. 

In general, surface interactions also perturb the degree of nematic ordering, e.g., the; value 
of S. Above the bulk nematic-isotropic transition temperature, this effect may induce the 
wetting of the surface by a nematic layer with S ^ 0, as predicted by Sheng [5]. This effect, 
long ago observed by Miyano [6], was also recently detected by a surface force apparatus 
measuring the related capillary condensation [7] . Although nematics are uniaxial in the bulk, 
the wetting layer may exhibit biaxiality due to the lower symmetry near the surface [8]. In 
the biaxial case, the symmetric and traceless Qij has its three eigenvalues {i = 1,2,3) all 
different. A nematic director, n, can still bo defined as the eigenvector corresponding to the 
eigenvalue of largest modulus, say A3. By analogy with the uniaxial case, one can then define 
the scalar ordcr-paramctcr by A3 = ^S; the amount of biaxiality, < B < 1, is then defined 
by Ai = -^^(l - B) and A2 = -^S{1 + B). 

The aim of the present paper is to model the planar degenerate anchoring of nematic 
liquid crystals, using the tcnsorial Landau-de Gennes formalism. To our knowledge, this has 
not been done before. Actually, there is some confusion in the literature, since a number of 
papers dealing with "degenerate planar anchoring" describe in fact homeotropic anchoring 
with negative S [8-11]. The confusion arises from the fact that the molecules themselves 
are planar and randomly oriented in this homeotropic situation. Again, conventional planar 
degenerate anchoring describes the tendency of a substrate to favor the common alignment of 
the molecules along any direction parallel to its surface [1] . This is the very situation that we 
consider in this paper. 

Theory. - Let us recall some of the existing tensorial anchoring models. The simplest 
quadratic surface free energy favoring a given non- degenerate nematic ordering Q^^^ i.e., a 
well-defined director and a well-defined scalar and, possibly, biaxial order-parameter, was 
proposed by Nobili and Durand [12]: 

/nd = \w (q,, - Qf) (q,, - Qf) . (1) 

Here, and throughout, summation over repeated indices is implicit. Because /nd is propor- 
tional to the sum of the squares of the components of Q — Q^'^\ it obviously has, for W > 0, 
a unique minimum for Q — Q'^^^ In this simple situation, the surface energy itself possesses 
a well-defined minimum. In principle, this is not necessary: the only requirement is that 
the nematic free energy as a whole, i.e., the sum of the surface and the bulk free energies, be 
bounded from below. For instance, early descriptions have dealt with only a linear surface free 
energy I'iQijVj, where v is the normal to the interface [5]. Actually, it is not clear at which 
order the surface free energy expansion should be truncated. In principle, for any truncation, 
one should verify that the higher-order terms are indeed negligible (this is especially critical in 
the case of large surface couplings). Another difficulty, when dealing with truncated surface 
free energies not bounded from below, is that one cannot determine which kind of anchoring 
they favor without including bulk terms and solving the whole problem. For these reasons, 
we shall concentrate in this paper on surface free energies that possess well-defined minima 
by themselves. 
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Fig. 1 - Polar contour plot of the surface potential 0, with Wi — W2 and ~ 0.5, for a uniaxial 
surface order-parameter Qij — S{ninj — |(5,;j), where n — (cos S cos 0, cos S sin <j!>, sin The angle 9 
is the inclination of the director with respect to the substrate and the radial coordinate corresponds 
to the scalar order-parameter — ^ < S < 1. The energy does not depend on (p. Darker shadings 
corresponds to lower energies. The absolute minimum is attained for degenerate planar anchoring: 
6 = (or tt) and S = So; the apparent homeotropic minimum for 9 = ir/2 and S < is actually a 
saddle-point with respect to the biaxial degree of freedom (outside the plane of this figure). 



The most general quadratic form describing isotropic substrates is [9, 13] 

/2 = Wii Vi Q.y Vj + 14^21 Qij Qij + W22 {vi Qij Vjf + 14^23 Qij Qjk Vk- (2) 

It contains all the scalars up to quadratic order made with Q and u (the surface normal). 
Minimizing /2 with respect to the five components of the most general symmetric traceless 
tensor yields, in the basis where the z-axis is along i^, a unique extremum [14], uniaxial, 
at Q(") = diag(A, A, -2A), where A = ^"11/(41^22 + 4W23 + 6W21). This corresponds to a 
homeotropic anchoring (n = u) with either 5 > or S" < depending on the sign of A. 

To obtain a planar degenerate uniaxial minimum, it is necessary to include terms of order 
higher than two. We shall go up to fourth order to have a surface free energy bounded from 
below. Note that this is also the lowest order for the bulk free energy to correctly describe 
the nematic-isotropic phase transition. Our aim is to obtain a degenerate manifold of minima 
for Qij = So {niTij — ^(^y), where n is an arbitrary unit vector parallel to the substrate, and 
where 5'o > is fixed. To this purpose, we define 

Qij = Qij + -^Sq Sij, (3) 

and we look for the constraints ensuring that Qij takes the form 5*0 rii rij , where n is parallel 
to the substrate. The two following conditions are necessary and sufficient [15]: i) Qij must 
coincide with its projection on the substrate Qj^ = Pik Qki Pij, where Pij = Sij — ViVj, and 
ii) the trace of the square of Qij must be equal to 5o^, that is Qij Qij — Sq'^. The simplest 
surface free energy being minimum if, and only if, these two conditions are satisfied is 

h = Wi (g,, - Qij) (q,, - Qij) + W2 (q.j Q^J - So^) ' , (4) 

with Wi > and W2 > 0. This surface potential is of fourth-order in Q and covariant. It 
is not, however, the most general fourth-order expansion giving rise to a planar degenerate 
anchoring: one can add, for instance, a contribution proportional to the square of the first 
term and several other similar, but linearly independent, terms. In the basis where the z-axis 
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is along the surface normal u, the above surface potential takes the explicit form 

2' 



/4 = Wi 2{Ql, + Ql^)+i^Q,, + Qyy 

+ 4 W2 (qI, + Qly + Qx. Qyy + + QL + 



5o 
3 



(5) 



where we have taken into account that Q is symmetric and traceless {Qzz — ~Qxx ~ Qyy)- 
As shown in fig. ^ the surface potential /4 possesses a degenerate absolute minimum for 
uniaxial order-parameters with planar director and 5* = 5*0 > 0. No other minima are present: 
the apparent homeotropic minimum with negative scalar order-parameter visible in fig. ^ is 
actually a saddle point with respect to biaxial perturbations. 

Wetting and anchoring properties. - To test the physical properties of the surface poten- 
tial Q), we consider a semi-infinite sample filled, for z > 0, with a liquid crystal either in the 
nematic or in the isotropic phase. For the bulk free energy density we take the fourth-order 
Landau-de Gennes expansion [1]: 



fv = (T - T*) Qij Qij 



1 1 2 

QijQjk Qki ~^ {Qij Qij^ 



2 Qij,k Qij,k : 



(6) 



where the comma indicates spatial derivation. The coefficients a, b, c are assumed positive, 
L > is an elastic constant (in the one constant approximation), T is the temperature, and T* 
is the supercooling temperature of the isotropic phase. We define the normalized temperature 
T = 27 acb~^ {T — T*) and the rescaled order-parameter qij = ^^/6b~^cQij . The first-order 
nematic-isotropic transition occurs for t = 1. In the nematic phase (r < 1), the reduced 

V9 — 8r). We rescale the lengths with respect 



scalar order-parameter is Sb(r) — 
to the correlation length at the transition ^ = b^^ \/ 27 cL and the energies with respect to Fq 
26[L/(3c)]'^/^. The dimensionless bulk and surface free energy densities become \t qij qij 



VGqi]qjk qki + ^{qtjqtj)'^ + ^qij,kqtj,k and wi{qij-qjj)[qij-qij)+W2{qij qij-so^Y, respectively, 
with the normahzed anchoring parameters wi = 36~^(3c/L)-^/^ Wi, W2 = |6(3Lc'^)~^/^ W2 
and So — ^V6b~^cSo- We look for equilibrium order-parameter profiles depending only on 
z and possessing a symmetry plane orthogonal to the substrate, which we take as the (y, z) 
plane. Therefore, we set 



q(^) 




(7) 



The minimization problem defined by fv and yields the Euler-Lagrange equations in the 
bulk: 



dz^ 
dz"^ 
dz^ 



V6( 



a 



2/32 + 7^ 



2a(3) + Aa (a^ 



T(3 



7 



76 

- 3%/6(c 



a/3), 

- 2a^ + 7^ - 2a(3) + 4/3 (a^ + ^32 _^ ^2 _^ a(3) , 
+ /3) + 4 + /32 + 72 + a(3) 



with the boundary conditions on the substrate, i.e., at z = 0, 
^ = ^W2a [3 (a^ + /32 + 72 + a/3) - sg] 



\wi (3/3- 



(8) 
(9) 
(10) 

(11) 
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Fig. 2 - a) Phase diagram of the surface state in the plane of the reduced temperature r and of 
the reduced anchoring strength w = wi = W2, for the preferred surface scalar order-parameter 
So = 1 < Sb{l). The states and H„ correspond to a planar biaxial ordering with positive scalar 
order-parameter and to a homeotropic uniaxial ordering with negative scalar order-parameter, re- 
spectively. Full (resp. dashed) lines indicate first (resp. second) order transitions. The second-order 
line terminates at a critical end-point, b) Same as a) but for so = 1.5 > Si,(l). The dot indicates 
a tricritical point. Inset: Evolution of the reduced surface scalar order-parameter, Ss, and of the 
amount of biaxiality, Bs, as a function of w, close to the second-order transition at r = 2. The latter, 
indicated by the dashed vertical line, actually occurs between H„ and a homeotropic negative biaxial 
state _ff^. At the maximum of biaxiality (dotted vertical line), the state transforms continuously 
into the state P\ the two phases having the same symmetry. 



dz 
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(13) 



and for z — > oo, da/dz — df3/dz — d'y/dz — 0. 

We solve numerically the above boundary value problem by a finite difference algorithm 
with deferred correction and Newton iteration [16]. When multiple solutions are present, we 
select the one with the lowest total free energy. For the sake of simplicity, we set wi = W2 = w 
and, for a given value of sq, we explore the behavior of the order-parameter profile for different 
values of r and w. We find that the eigenvectors of the order-parameter solutions are always 
parallel to the a;, y and z axes, which corresponds to 7(2:) = 0. 

In the nematic phase (r < 1), we find the expected planar anchoring, with, however, a 
slight biaxiality in a surface layer of thickness ~ ^. The nematic director n is parallel to 
the surface and uniform along z, the scalar order-parameter s > relaxes over ~ ^ to the 
bulk value Sf, starting from a surface value Ss intermediate between sq and Sb- The amount 
of biaxiality at the surface, Bs, vanishes for sq = Sb and grows with |so ~ Sfeh it obviously 
vanishes for w — 0, but also for w —^ 00, presenting a maximum at a finite value of w. 

In the isotropic phase (r > 1), the order-parameter vanishes in the bulk, but a wetting 
nematic layer appears at the surface. For large values of w, the ordering at the surface is 
planar and biaxial, as in the nematic phase. We call this state (the dagger indicating 
biaxiality). As before, the amount of biaxiality tends to zero as w — > 00. Conversely, for small 
values of w, the ordering at the surface is homeotropic and uniaxial, with a negative scalar 
order-parameter. We call this state _ff„. In between, a surface transition occurs. 

In fig. 121 we show the typical surface phase diagrams in the plane (r, w) for two different 
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values of sq. At the bulk phase transition r — 1, for small values of w there is always a 
first-order transition between the states and iJ„. For larger values of w, two possible 
scenarios occur, depending on the value of sq with respect to the bulk value at the transition 
Sb(l) = •\/3/2 ~ 1.22. When sq < s;,(l) (fig- there is a first-order transition between 
two different states. When sq > Sb(l) (fig-Eb), no surface transition occurs: approaching 
T = 1 from above, a wetting nematic layer of diverging thickness continuously invades the 
whole sample. This corresponds to a complete wetting situation, similar to the uniaxial case 
described in Ref. [5,6,17]; indeed, the surface-induced biaxiality relaxes over ~ ^. In the 
isotropic phase (t > 1), a surface phase transition occurs at a given value of w. For sq > st 
and T < Tc (corresponding to the tricritical point in fig. EId), the transition is of first order, 
between the states and . For sq < Sb (fig-Et); and for t > Tc when sq > si, (fig.Eb), the 
transition is of second order, between i7„ and a biaxial homeotropic negative state, H^. In 
fact, by symmetry, it is impossible to go continuously from _ff„ to P^ , since, when the nematic 
director switches from one eigenvector to another, the amount of biaxiality must be maximum 
[the order-parameter tensor is then of the form diag(0. A, — A)]. This behavior is illustrated in 
the inset of fig. Eb- Note that the transformation from H}^ to P^ is a continuous evolution 
between two states having the same symmetry and not a phase transition: the discontinuity 
in Ss and the cusp in B^, which occur when the biaxiality is maximum, do not correspond 
to any actual discontinuity of q. In fact, close to this transformation, the tensor q{z) at the 
surface z = can be written as 

/e \ 
q(0) =\0 a-e , (14) 
\0 -a J 

with cr > and e a small quantity (|e| < a) continuosly varying from positive to negative 
across the transformation. In the state, e > 0: the eigenvalue of q(0) of maximum 
modulus is —a < 0. This corresponds to a nematic director along z (homeotropic ordering), 
with the negative scalar order-parameter Ss = —3a/2 and the amount of biaxiality Bs = 
1 — 2e/a. In the state, e < 0: the eigenvalue of maximum modulus is now cr — e > 0. 
This corresponds to a nematic director along y (planar ordering), with the positive scalar 
order parameter Ss — 3(cr — e)/2 and the amount of biaxiality Bg — 1 + 2e/{a — e). At the 
transformation, e = 0: no nematic director can be defined and the amount of biaxiality is 
maximum (Bs — 1). Therefore, the discontinuous change of sign of and the cusp of Bs at 
the transformation are not due to a discontinuity of q(0) but to the change of the direction 
associated to the eigenvalue of maximum modulus. Note also that the width of the region 
is extremely narrow and goes to zero at the tricritical point. 

Our findings are compatible with the experimental data of Tarczon and Miyano [17], who 
observed a pretransitional uniaxial negative birefringence with optical axis perpendicular to 
the surface, corresponding to our Hn state, for the compound MBBA in contact with a silane- 
treated substrate giving planar degenerate anchoring in the nematic phase. This negative 
pretransitional birefringence was not observed for the compound 5CB on the same substrate, 
which probably corresponds to our P^ surface state. 

Conclusions. - As we have shown, the surface free energy describing a planar degenerate 
anchoring in the framework of the Landau-de Gennes formalism must be at least of fourth- 
order in the tensorial order-parameter Q. Here, we have proposed the simplest expression 
having this property and we have investigated the associated anchoring and wetting behavior. 
In the nematic phase, we find the expected planar degenerate anchoring with, however, a small 
surface biaxiality. In the isotropic phase, a surface wetting layer appears: for large enough 
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surface couplings it is planar biaxial, while for low couplings it is homeotropic uniaxial with 
negative scalar order-parameter. These predictions are compatible with the experimental 
observations for substrates yielding planar degenerate anchoring [17]. The transition between 
the two surface states may be of first or second order. The biaxiality in the planar states 
tends to zero for large surface coupling, yielding a uniaxial degenerate planar anchoring or 
wetting. It would be interesting to study the consequences of this degenerate anchoring in 
more complicated geometries such as, e.g., colloidal suspensions. 

Note added in proofs. - After this work was completed. Prof. Sluckin kindly pointed 
out to us that the surface potentialof ref. [8], linear in the tensor order-parameter, gives a 
genuine planar degenerate anchoring in the case of complete wetting. We checked numerically 
that this is true, although the ordering is strongly biaxial in a coherence length close to the 
surface. Ac;cordingly, above the transition temperature, where the wetting layer is of thickness 
comparable to the coherence length, the surface layer is strongly biaxial. Our fourth-order 
potential, on the other hand, gives an almost perfectly planar ordering in the whole microscopic 
layer, since it directly favors the planar anchoring. 

* * * 

We thank B. Silva for a critical reading of the manuscript and useful comments. 
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